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■  Investigates  the  design  of  a  transpolar  extremely-low-frequencv  (ELF) 
experiment  that  would  validate  Defense  Nuclear  Agency  propagation  codes  for 
ELF  system  performance  in  nuclear  environments.  To  he  conducted  during  the 
next  solar  maximum,  when  polar  cap  absorption  events  (PCAs)  might  approximate 
a  nuclear  environment,  the  ideal  experiment  would  simultaneously  measure  both 
propagation  and  ionospheric  structure.  Propagation  measurements  would  be  taken 
at  receivers  located  such  that  the  direct  propagation  path  did  not  approach 
the  polar  cap  boundary;  the  ELF  values  would  then  be  compared  with  computer 
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» code  outputs.  A  moderate  PCA  might  cause  a  3  to  4  dB  reduction  in  a  transpolar 
ELF  signal;  a  strong  event,  an  8  or  9  dB  reduction.  At  least  2  hours  would  be 
required  to  resolve  amplitude  changes  in  the  signal;  noise  processing  to  excise 
large  pulses  would  be  essential.  Measurements  of  ionospheric  changes  would  be 
taken  with  rocket-borne  instruments,  entered  into  computer  calculations  using 
the  DNA  codes,  then  compared  with  the  propagation  measurements.  Rocket  data 
should  be  gathered  at  50  km  or  below;  data  from  previous  PCAs  obtained  only 
at  60  km  or  above  yield  unacceptable  uncertainties  in  calculated  field  strengths 
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SUMMARY  AND  CONCLUSIONS 


This  report  investigates  the  design  of  a  transpolar  extremely- 
low-frequency  (ELF)  experiment  that  would  validate  Defense  Nuclear 
Agency  codes  for  predicting  ELF  propagation  in  nuclear  environments. 
Such  an  experiment  would  be  conducted  during  a  maximum  in  the  solar 
cycle,  when  disturbances  in  the  polar  ionosphere — polar  cap  absorp¬ 
tion  events,  or  PCAs — caused  by  solar  disruptions  would  resemble 
distortions  in  the  earth-ionosphere  waveguide  caused  by  radiation 
from  a  high-altitude  nuclear  burst. 

The  most  effective  experiment  would  simultaneously  measure  both 
propagation  and  ionospheric  structure.  Propagation  measurements 
would  be  accomplished  with  one  or  more  receivers  placed  (e.g.,  in 
Norway)  so  that  an  ELF  signal  transmitted  from  the  Wisconsin  Test 
Facility  (WTF)  traversed  at  least  a  few  megameters  of  the  polar  cap 
(but  without  obliquely  intersecting  the  cap  boundary,  to  avoid  trans¬ 
verse  edge  effects  that  would  complicate  interpretation  of  "the  measure¬ 
ments).  Changes  in  the  WTF  signal  measured  during  a  PCA  would  then 
be  compared  with  field-strength  values  calculated  from  the  ELF  pre¬ 
dictive  codes — thereby  checking  the  codes.  Data  necessary  for  cal¬ 
culating  refractive-index  height  profiles  (the  basic  inputs  to  the 
codes) — height  profiles  of  electron  density,  ion  density,  or  positive 
and  negative  conductivity  in  the  lower  ionosphere — would  be  obtained 
using  rocket-borne  instruments 

The  propagation  measurements  would  have  to  determine  signal 
changes  whose  uncertainty  was  much  smaller  than  the  expected  PCA- 
induced  reductions  (3  to  4  dB  from  a  moderate  PCA,  8  or  9  dB  from  a 
strong  PCA);  otherwise,  experimental  errors  would  mask  the  PCA's 
effect.  Resolving  amplitude  changes  from  a  large  PCA  would  require 
at  least  two  hours  of  integration,  since  a  transpolar  WTF  signal 
would  be  very  weak.  Noise  processing  to  excise  large  pulses  would 
be  essential. 

To  minimize  uncertainties  in  ELF  signal  strength  calculated  from 
the  predictive  codes,  electron  density  data  should  be  obtained  at 
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altitudes  of  50  km  or  below,  where  electrons  strongly  affect  ELF 
propagation;  rocket  measurements  of  ionospheric  structure  made  during 
previous  PCAs  obtained  data  only  for  altitudes  of  60  km  or  above, 
yielding  unacceptable  uncertainties.  To  obtain  data  at  lower  alti¬ 
tudes  would  require  sensing  a  few  tens  of  electrons  per  cubic  centi¬ 
meter  in  a  collision-dominated  region. 

Uncertainties  in  interpreting  the  results  (that  is,  in  comparing 
measured  with  calculated  changes  in  field  strength)  related  to  the 
excitation  factor  at  the  ELF  receiver  could  be  avoided  by  (1)  locating 
the  receiver  so  that  there  was  no  doubt  whether  it  was  inside  or  out¬ 
side  the  irradiated  cap,  or  (2)  sounding  the  ionosphere  at  the  re¬ 
ceiver  site. 
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PREFACE 


This  report  documents  contingency  planning — requested  by  the 
Defense  Nuclear  Agency — for  a  transpolar  extremely- low- frequency 
(ELF)  propagation  experiment  to  be  performed  during  the  forthcoming 
maximum  in  the  solar  cycle.  The  experiment  would  verify  theoretical 
predictive  codes  under  disturbed  ionospheric  conditions.  Even  though 
budgetary  constraints  will  probably  prohibit  the  experiment,  the 
present  analysis  should  assist  the  designers  of  future  experiments 
that  require  coordinated  measurements  of  ELF  signal  strength  and 
ionospheric  structure. 
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SECTION  1 
INTRODUCTION 

Propagation  codes  developed  by  the  Defense  Nuclear  Agency  (DNA) 
to  predict  extremelv-low-f requency  (ELF)  system  performance  in  nuclear 
environments  have  not  been  fully  verified  experimentally.  Pacific- 
Sierra  Research  Corporation  (PSR)  was  requested  by  DNA  to  investigate 
certain  aspects  of  the  design  of  a  transpolar  propagation  experiment 
to  validate  these  codes.  Such  an  experiment  would  be  performed  during 
a  maximum  in  the  solar  cycle,  when  polar  cap  absorption  events  (PCAs) 
would  be  expected  to  simulate  salient  aspects  of  nuclear  environments. 
This  report  summarizes  the  results  of  our  investigation  and  establishes 
minimum  tolerances  for  a  definitive  verification  of  ELF  propagation 
codes . 

Radiation  from  high-altitude  nuclear  bursts  ionizes  the  atmosphere 
below  the  ambient,  ionospheric  D-  and  E-layers — a  distortion  of  the 
earth-ionosphere  waveguide  that  may  affect  military  long-wave  communica¬ 
tion  systems.  Computer  codes  predict  that  a  link  with  at  least  several 
megameters  of  propagation  path  exposed  to  such  an  environment  can 
suffer  significant  signal  loss.  That  conclusion  has  been  substantiated 
experimentally  for  low-  and  very-low-frequencv  (LF/VLF)  links,  for 
which  data  are  available  from  (1)  the  Fish  Bowl  high-altitude  nuclear 
test  series,  (2)  a  laboratory  model  of  the  earth-ionosphere  waveguide 
[Field  et  al . ,  1972],  and  (3)  propagation  experiments  under  naturally 
disturbed  conditions  [Westerlund  et  al.,  1969]. 

ELF  communications  (45  to  100  Hz)  are  much  newer  than  LF/VLF  com¬ 
munications,  operational  for  decades.  The  only  existing  ELF  transmitter 
is  the  low-power  Wisconsin  Test  Facility  (WTF)  constructed  quite 
recently.  The  data  on  ELF  propagation  are  thus  much  less  comprehen¬ 
sive  than  those  on  LF/VLF.  Indeed,  data  on  disturbed  conditions  are 
limited  to  a  few  fortuitous  measurements  of  natural  ELF  signals  (atmo¬ 
spheric  noise)  during  high-altitude  tests  and  natural  ionospheric 
d i sturbances . 
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Field  [1978]  found  the  sparse  data  generally  consistent  with  com¬ 
puter  predictions  of  degraded  propagation  for  severe  and  widespread 
disturbances.  Nonetheless,  no  experiment  has  yet  monitored  a  signal 
from  a  controlled  ELF  source  during  a  strong  ionospheric  disturbance. 
The  probable  increase  in  the  number  and  magnitude  of  P(  ,s  curing  the 
approaching  maximum  in  the  11  year  solar  cycle  should  afford  such  an 
opportunity . 

A  PCA  approximates  certain  features  of  a  nuclear  environment, 
since  solar  conditions  cause  energetic  protons  to  penetrate  and  ionize 
the  polar  ionosphere  to  altitudes  as  low  as  40  km.  The  most  effective 
experiment  would  simultaneously  measure  propagation  and  ionospheric 
structure.  Propagation  measurements  would  quantify  changes  in  the  K'TF 
signal  during  a  PCA,  thus  providing  values  to  compare  with  computer  cod 
outputs.  Such  measurements  would  be  accomplished  with  one  or  more  re¬ 
ceivers  placed  so  the  signal  traverses  at  least  a  few  megameters  of 
the  polar  cap.  Values  quantifying  changes  in  the  ionospheric  struc¬ 
ture  would  be  entered  into  the  computer  programs  for  calculating  field- 
strength  changes  to  compare  with  measured  values — thereby  checking 
the  codes.  Such  measurements  could  be  accomplished  with,  for  example, 
rocket-borne  instruments. 

Our  analysis  shows  that  a  moderate  PCA  might  cause  a  3  to  4  dB  re¬ 
duction  in  a  transpolar  ELF  signal,  whereas  a  strong  event  might  cause 
an  8  or  9  dB  reduction.  The  propagation  measurements  must  determine 
signal  changes  with  an  uncertainty  much  smaller  than  the  expected  PCA- 
induced  reductions;  otherwise,  experimental  errors  will  mask  the  effect 
to  be  studied.  Similarly,  ionospheric  measurements  used  as  code  inputs 
must  be  accurate  enough  for  meaningful  comparison  of  calculated  and 
measured  signal  changes. 

Sections  2  and  3  consider  the  measurement  of  ionospheric  struc¬ 
ture.  Section  2  examines  data  from  rocket  samplings  of  vertical  iono¬ 
spheric  structure  during  previous  PCAs,  identifying  inadequacies  tiiat 
must  be  corrected  for  future  experiments.  Section  3  examines  the 
propagation  effects  of  lateral  ionospheric  structure  in  the  polar  cap, 
quantifies  computational  uncertainties  caused  by  such  structure,  and 
compares  the  uncertainties  with  the  minimum  tolerances  required  for 
code  verification. 


Section  4  addresses  the  problem  of  weak  transpolar  WTF  signal, 
which  must  be  integrated  for  an  accurate  estimate  of  received  signal. 
Too  short  an  integration  time  causes  unacceptable  uncertainty  in 
amplitude.  Conversely,  too  long  an  integration  time  may  create  un¬ 
certainties  due  to  substantial  changes  in  the  incident  flux  and  iono¬ 
spheric  ionization  during  measurement.  Criteria  are  established  for 
the  best  compromise  on  integration  time.  The  benefits  of  nonlinear 
suppression  of  atmospheric  noise  at  the  receiver  are  also  quantified. 
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SECTION  2 


ROCKET  MEASUREMENTS  OF  IONOSPHERIC  STRUCTURE  FOR  USE  AS 
INPUTS  TO  ELF  PROPAGATION  CALCULATIONS 

Rocke "-borne  instruments  measure  the  height  profiles  of  electron 
density,  ion  density,  or  positive  and  negative  conductivity  in  the 
lower  ionosphere.  Such  data  are  needed  to  calculate  refractive-index 
height  profiles,  which  are  the  basic  inputs  to  ELF  predictive  codes. 
Thus  rocket  data  permit  a  direct  comparison  of  measured  and  calculated 
ELF  field  strengths  for  ionospheric  conditions  at  the  time  of  propaga¬ 
tion.  Here  we  summarize  the  formulas  for  converting  measured  iono¬ 
spheric  parameters  to  refractive  index  profiles,  identify  inaccuracies, 
and  establish  criteria  for  precision  of  measurement. 

FORMULAS  FOR  REFRACTIVE  INDEX 

Rocket  data  can  most  easily  be  converted  to  the  refractive  index 
when  the  positive,  and  negative  conductivities  0+  and  a_  are  measured 
directly.  In  this  case,  the  refractive  index  n  is  calculated  from 

n2  =  1  -  -f-  (O+O)  .  (1) 

weo  + 

Since  =  8.85  x  10  ^2  and  the  angular  frequency  u>  is  known  exactly, 
the  accuracy  of  the  refractive  index  computed  from  Eq .  (1)  is  equiva¬ 
lent  to  that  of  the  measured  data  for  the  conductivities.  However, 

Eq .  (1)  is  valid  only  if  <J+  and  O  denote  the  cumulative  conductivities 
of  all  positive  and  negative  species,  respectively.  It  could  he 
grossly  in  error  at  certain  altitudes  if,  for  example,  the  instrument 
used  to  measure  a  could  fully  sense  the  conductivity  of  ions  hut  not 
electrons. 

Unfortunately,  the  refractive  index  is  less  accurate  when  based 
on  electron  density  N,  positive  ion  density  N+,  and  negative  ion 
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density  instead  of  conductivity.  In  that  situation,  the  ELF  re¬ 
fractive  index  at  heights  below  about  75  km  is  calculated  from 


(2) 


where  q  is  the  electron  charge,  Mg  is  the  electron  mass,  V  is  the 
electron  collision  frequency,  is  the  mass  of  the  ath  positive  ion, 
V  the  collision  frequency  of  the  ath  positive  ion,  and  so  forth. 
Comparing  Eq .  (2)  with  Eq .  (1),  we  find  that 


a 


+ 


and 


a 


Neq2  y'  V2 

M  V  +  Z-f  Mqvq  ‘ 
e  e  6  6 


(3a) 


(3b) 


Tlius  the  refractive  index — or,  equivalently,  the  conductivities — must 
be  synthesized  by  summing  over  all  charged  species  sufficiently  popu¬ 
lous  to  make  a  significant  contribution.  The  terms  in  the  sums  are 
not  governed  solely  by  the  particle  densities,  but  instead  depend  on 
the  combination  N/Mv  for  each  species.  That  dependence  causes  no 
serious  loss  of  accuracy  when  the  conductivity  (or  refractive  index) 
is  dominated  bv  electrons,  because  the  electron  mass  is  exact  and  V 

e 

is  believed  to  be  known  accurately.  However,  the  ion  collision  fre¬ 
quencies  are  uncertain  by  a  factor  of  2,  and  masses  of  important  ion 
species  are  often  imprecise.  At  lower  altitudes  where  ions  dominate 
the  refractive  index,  therefore,  thi’  aombi tuition  N/Mv  if  known  no 


To  accommodate  the  large  number  of  ionic  species,  we  use  the 
subscripts  a  and  B  to  denote  the  "ath"  and  "flth"  positive  and  nega¬ 
tive  species. 
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u>cter  than  to  within  a  factor  of  2  or  3  even  if  the  ion  density  N 
is  precise.  In  principle,  that  uncertainty  can  be  avoided  by  directly 
measuring  conductivities  at  ion-dominated  altitudes,  because,  as  Eq . 

(3)  shows,  such  measurements  directly  determine  the  sum  of  N/MV  over 
all  important  species. 

In  practice,  the  full  sums  in  Eq.  (2)  are  seldom  used  in  engineer¬ 
ing  calculations.  Ionic  species  are  classified  as  simply  "positive" 
or  "negative"  ions,  each  type  having  properties  determined  by  averag¬ 
ing  over  the  species  actually  present.  Further,  it  is  usually  assumed 
that  such  generic  ions  have  equal  masses  and  collision  frequencies 

V^.  Equation  (2)  then  takes  the  simplified  form 


2 

n 


1  - 


N+  + 


M.V. 

l  i 


(A) 


where  N+  and  N  are  the  densities  of  the  "average"  positive  and  nega¬ 
tive  ions,  respectively.  The  values  of  M .  and  v.  at  a  given  alti¬ 
tude  are  uncertain.  We  use  an  atomic  weight  of  32  here  for  M., 
but — depending  on  altitude — values  as  low  as  19  or  as  high  as  50  or 
60  are  possible.  Similarly,  V  could  be  anywhere  from  a  tenth  to  a 

fortieth  of  V  ;  we  follow  current  convention  in  assuming  V.  =  V  /20. 

e  i  e 

We  use  Eq.  (1)  to  compute  the  refractive  index  from  blunt-probe 
conductivity  data.  Inserting  the  value  for  £q  and  assuming  a  fre¬ 
quency  of  75  Hz,  Eq .  (1)  becomes 


n2  -  1  =  -2 . 4i  x  io8  (cr+  +  0_)  ,  (5) 

where  and  0  are  in  mhos  per  meter. 

We  use  Eq.  (4)  to  compute  the  refractive  index  from  data  on 
electron  and  ion  density.  Since  N  R*  N+  whenever  ions  contribute 
significantly  to  the  refractive  index,  we  can  rewrite  Eq .  (4)  as 


2 

n 


1  = 


(!)£_.  M  V 
0  e  e 


f  2M  V  A 

N  +  N  ) 

e  M.V.  I 

\  1  1  / 


(6) 


* 

6 
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Substituting  numerical  values  and  assuming  =  32  and  v.  =  v^/20, 
F.q.  (6)  becomes 


2 

n 


1 


(7) 


where  and  N  are  in  units  of  (meters)  and  is  taken  from  Knapp 

and  Schwartz  [1975,  Fig.  5-1]. 

2 

We  can  determine  n  accurately  from  density  data  at  high,  electron- 
dominated  altitudes  because  the  parameters  in  the  electron  contribu¬ 
tions  to  Eq.  (7)  are  well  known.  Similarly,  at  low,  ion-dominated 
altitudes,  the  problem  of  poorly  known  and  NT  can  be  circumvented 
by  measuring  ion  conductivity  directly  and  using  Eq .  (5)  to  determine 
the  refractive  index.  Potential  difficulties  arise  in  the  transition 
region  where,  although  much  less  numerous  than  ions  and  difficult  to 
measure,  electrons  still  make  a  major  contribution  to  n^.  Unfortun¬ 
ately,  that  region  occurs  at  altitudes  crucial  to  EI.F  propagation. 

For  example,  as  Eq.  (7)  demonstrates,  the  transition  from  electron 
to  ion  dominance,  occurs  at  the  altitude  where 


n  *Jii- 

e  1500  ’ 

More  explicitly,  we  can  write 

+ 

N~ 

N ^  >,  2~5q~q  (electrons  dominate) 

and 

+ 

N" 

Ne  ~  1500  (i°nS  dominat°)  ■ 


(Ha) 


(8b) 


(He) 


Even  at  altitudes  where,  say,  ~  10  N  ,  therefore,  electron  density 
must  be  measured  accurately.  As  shown  below,  this  requirement  im¬ 
plies  measurement  of  a  few  tens  of  electrons  per  cubic  centimeter  at 
altitudes  of  50  to  60  km  during  a  PCA. 


r< 


r.. 


!, 
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EXAMPLE:  FCA  09 

Using  data  from  the  PCA  of  2  November  1969,  we  now  consider 
problems  encountered  in  using  rocket  measurements  of  ionospheric 
structure  for  ELF  propagation  calculations.  Other  PCAs  have  been 
stronger  or  possibly  more  interesting,  but  none  has  been  so  carefully 
monitored.  A  series  of  rocket  measurements  using  several  sensor  types 
are  reported  in  the  Ptwcedi  njp  of  ( '< :!'A  H  i'lit’ij-oriuir  im  I'uvr'olo 

Fi'rnt  of  NoVi,rnbrr‘  1909,  henceforth  called  !'<'A  09.  The  papers  bv  Hale 
et  al.  and  Ulwick,  cited  below,  appear  in  that  report. 

Of  primary  interest  are  blunt-probe  measurements  of  o  and  o 
given  by  Hale  et  al.  [1972]  and  the  following  measurements  given  by  1 

Ulwick  [1972]:  Langmuir  probe  (LP)  measurements  of  N>  and  N+,  Z-0 
probe  measurements  of  ,  and  Gerdien  condenser  (GC)  measurements  of 
N+.  Here  we  compute  the  ELF  refractive  index  using  selected  samples 
of  those  data,  thereby  stressing  inconsistencies  among  the  measure¬ 
ments  and  altitude  gaps  where  the  data  are  inadequate. 

2 

Table  1  lists  |n  -  l|  at  various  altitudes,  calculated  from 
Eq .  (5)  using  the  conductivity  measurements  of  Hale  et  al.  [1972]; 
and  from  Eq .  (7)  using  the  electron  and  ion  density  data  summarized 
by  Ulwick  [1972].  Results  are  given  for  sunrise,  daytime,  and  sunset 
on  3  November  1969 — the  day  after  the  onset  of  the  event. 

For  convenience.  Figs.  1  through  3  plot  the  refractive  index 
2 

curves  for  |n  -  l|  at  various  altitudes  from  Table  1.  The  lower 
altitude  portions  of  the  curves  show  ions  only,  since  Ulwick  reports 
no  N  values  for  those  regions  (recall  that  N  +  N  ~  2N  )  .  At 
sufficiently  low  altitudes,  the  "ions  only"  refractive  index  should 
approximate  the  total  index.  At  higher  altitudes,  the  Ulwick  data 
give  the  total  refractive  index,  including  botli  electrons  and  ions, 
but  are  dominated  by  the  highly  mobile  electrons.  At  altitudes  for 
which  Table  1  lists  data  for  both  N+  and  N  ,  we  have  prorated  t  lie 
contribution  of  ions  according  to  Eq .  (7),  thereby  extending  the 
"ions  only"  curves. 
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Table  1.  Comparison  of  refractive-index  height  profiles  at  75  Hz 
computed  from  various  rocket  data;  3  November  1969. 
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NOTF:  Not  .ill  sensor*  were  operating  at  the  three  time**  considered:  Male  «•  t  al. 

give  no  conduct  ivl  tv  data  for  sunset,  and  Ulwick  obtained  no  cerdirn  tondenser  mea¬ 
surements  at  sunrise.  All  measurements  reported  bv  Hale  et  al.  were  made  with  a 
blunt  probe,  whereas  the  method  (Z-H,  l.P,  or  Cl  )  used  to  obtain  the  data  given  bv 
Ulwick  are  indicated  at  the  top  of  each  column.  IT  »  local  time.  halt'  et  al  .  give 
values  of  conduct  lvi  tv  in  mhos  per  meter.  Ulwick's  measurements  are  in  unit*-  of 
(cent i meters ) " K 
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Figure  2.  Refractive-index  height  profiles  at  75  Hz:  daytime,  3  November  1969 


Refractive-index  height  profiles  at  75  Hz:  sunset,  3  November  1969 
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Several  conclusions  can  he  drawn  from  Table  1,  or  the  equivalent 
Figs.  1  through  3: 

1.  Above  60  km  or  so,  electrons  dominate  propagation  because 
the  total  refractive  index  calculated  from  the  conductivity 
data  is  much  larger  than  that  calculated  from  ions  only. 

2.  At  altitudes  below  about  55  km,  the  refractive  index  from 
Ulwick's  data  is  somewhat  smaller  than  that  from  the  con¬ 
ductivity  data.  Tiie  difference  could  simply  be  due  to 
scatter  between  two  measurement  methods,  or  could  occur  be¬ 
cause  the  conductivity  data  include  some  electron  contribu¬ 
tion  (as  Hale  et  al.  assert)  whereas  Ulwick's  data  do  not. 

3.  Large  gaps  in  the  refractive  index  profiles  occur  at  alti¬ 
tudes  between  about  60  km,  above  which  electron  density  data 
are  reliable;  and  45  to  50  km,  below  which  ions  make  the 
dominant  contribution. 

The  third  item  above — the  gaps  in  the  refractive  index  profiles 
between  about  50  and  60  km — probably  causes  the  greatest  uncertainty 
in  interpreting  ELF  propagation  experiments.  Such  gaps  occur  because 
electrons  contribute  strongly  down  to  the  altitude  where  ~  N  /1500. 
Since  ion  densities  of  a  few  tens  of  thousands  per  cubic  centimeter 
are  typical  for  PCA  conditions  (see  Table  1),  electron  densities  as 
small  as  a  few  tens  of  electrons  per  cubic  centimeter  must  be  mea¬ 
sured  for  a  good  estimate  of  the  refractive  index.  The  data  give 
electron  densities  no  smaller  than  several  hundred  per  cubic  centimeter — 
concentrations  that  occur  near  60  to  65  km.  Electron  data  are  needed 
for  altitudes  perhaps  10  km  lower,  where — though  much  less  numerous  than 
ions — electrons  are  nevertheless  the  dominant  influence. 

Therefore,  two  types  of  uncertainties  would  have  arisen  if  a 
coordinated  rocket/ELF-propagat ion  experiment  had  been  performed  dur¬ 
ing  PCA  (P).  First,  differences  exist  between  refractive  indexes  ob¬ 
tained  by  alternative  methods  in  the  ion-dominated  region,  which 
includes  most  altitudes  below  ~50  km.  Second,  the  interpolation  of 
values  in  the  refractive  index  between  the  ion-dominated  region  and  the 
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lowest  electron-dominated  altitudes  for  which  electron  data  are  avail¬ 
able  is  poorly  defined. 

To  quantify  these  uncertainties,  we  have  calculated  ELF  attenua¬ 
tion  rates  for  the  four  refractive  index  profiles  shown  in  Fig.  4. 

All  profiles  pertain  to  sunrise,  3  November  1969,  and  represent  four 
models  obtained  by  reasonable  fits  to  the  data  plotted  in  Fig.  1. 

They  are  distinguished  as  follows: 

•  Profile  A  represents  a  sharp  roll-off  from  the  index  based 
on  Ulwick's  data  in  the  electron-dominated  region  to  that 
based  on  his  data  in  the  lower  altitude,  ion-dominated  region. 

•  Profile  B  represents  a  transition  from  the  index  based  on 
Ulwick's  data  in  the  electron-dominated  region  to  that  based 
on  the  conductivity  data  of  Hale  et  al.  in  the  ion-dominated 
region. 

•  Profile  C  represents  a  gentle  roll-off  from  the  index  based 
on  Ulwick's  data  in  the  electron-dominated  region  to  that 
based  on  his  data  in  the  ion-dominated  region. 

•  Profile  D  is  analogous  to  profile  B,  but  was  chosen  to  have 
a  greater  height  gradient  between  altitudes  of  40  to  50  km, 
which  strongly  affect  ELF  propagation. 

Differences  among  the  profiles  are  due  to  the  inadequacies  of  the 
rocket  data  discussed  above. 

Table  2  lists  attenuation  rates  in  decibels  per  megameter  of 
propagation  calculated  with  PSR's  long-wave  propagation  code  using 
the  profiles  of  Fig.  4.  It  also  catalogs  the  calculated  PCA- induced 
losses  for  the  profiles,  based  on  the  assumptions  that  a  5  Mm  length 
of  propagation  path  is  exposed  to  the  i’CA  and  that  the  ambient  attenu¬ 
ation  rate  is  1.15  dB/Mm.  A  transpolar  propagation  experiment  would 
be  designed  to  measure  such  losses,  which  are  the  disparities  between 
the  received  signals  under  ambient  and  PEA  conditions.  The  differences 
among  the  computed  losses  result  from  both  scatter  in  the  rocket  data 
and  ambiguities  in  interpolating  the  refractive  index  through  altitudes 
ranging  from  about  45  to  60  km. 


Table  2.  Calculated  attenuation  rates  and 
PCA-induced  signal  losses  for 
75  Hz  and  5  Mm  of  exposed  path. 


Profile 

Attenuation 

Rate 

(dB/Mm) 

Signal 
Loss  (dB) 

A 

1.5 

1.75 

B 

1.8 

3.25 

C 

1.35 

1.0 

D 

1.5 

1.75 

For  a  useful  comparison  with  experimental  data,  the  computed 
results  should  exhibit  a  scatter  significantly  smaller  than  that  of 
the  PCA-induced  loss.  Unfortunately,  Table  2  shows  that  the  scatter 
is  a  large  fraction  of  the  expected  signal  loss.  Strictly,  this  con¬ 
clusion  applies  to  only  the  PCA  P9  event.  However,  similar  difficulties 
could  be  expected  for  any  event  where  the  accuracy  and  sensitivity  of 
rocket-gathered  measurements  are  no  better  than  for  PCA  t'i>.  Future 
experiments  should  concentrate  on  obtaining  electron  data  down  to 
altitudes  of  50  km  or  even  lower. 

Should  electron  data  prove  unobtainable  at  such  low  altitudes 
with  rocket-borne  instruments,  an  alternative  approach  would  be  needed. 
One  possibility  is — applying  atmospheric-chemical  relations — to  theo¬ 
retically  extrapolate  measured  electron  densities  to  altitudes  below 
60  km.  Another  is  to  use  satellite-borne  sensors  to  measure  fluxes 
of  ionizing  radiation  incident  on  the  ionosphere;  then,  using  the 
chemical  relations,  to  calculate  the  electron  density  at  altitudes 
ranging  from  45  to  60  km.  Both  alternatives  are  subject  to  uncer¬ 
tainties  in  the  values  of  ionospheric  reaction-rate  coefficients. 
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SECTION  3 


EFFECTS  OF  LATERAL  IONOSPHERIC  STRUCTURE  ON  ELF  DURING  PCA 

The  previous  section  considered  the  effects  of  vertical  iono¬ 
spheric  structure  on  ELF  propagation  during  a  PCA.  Here  we  consider 
lateral  ionospheric  structure,  which  can  influence  ELF  propagation  for 
two  reasons.  First,  unlike  conventional  frequencies,  ELF  cannot  be 
analyzed  solely  in  terms  of  ionospheric  properties  in  the  immediate 
vicinity  of  the  direct  propagation  path.  Instead,  ionospheric  proper¬ 
ties  must  be  averaged  over  an  ELF  Fresnel  zone,  which  can  cover 
millions  of  square  kilometers  for  a  transpolar  path.  Second,  the 
finite  size  of  the  polar  cap  must  be  taken  into  account. 

SATELLITE  DATA 

Ideally,  data  on  ionospheric  height  profiles  at  multiple  dis¬ 
persed  locations  are  needed  to  determine  the  extent  and  uniformity  of 
a  PCA- induced  disturbance.  To  our  knowledge,  unfortunately,  no  co¬ 
ordinated  multiple-location  measurements  of  ionospheric  conductivity 
or  ionization  have  been  taken  during  a  major  PCA.  Even  the  extensive 
rocket  measurements  during  PCA  69  were  taken  at  one  location.  Thus, 
data  on  the  lateral  structure  of  the  earth-ionospheric,  waveguide 
during  a  PCA  must  be  gleaned  from  satellite  measurements  of  incident 
particle  fluxes,  which  enhance  ionization.  Such  a  use  of  fluxes  in¬ 
stead  of  conductivity  or  density  profiles  introduces  uncertainties 
associated  with  the  chemical  calculations  used  to  obtain  the  ioniza¬ 
tion  levels. 

Reagan  and  Watt  [1976]  and  Reagan  et  al .  [1978]  give  data  on 
particle  fluxes  as  functions  of  latitude  for  the  PCAs  of  3  August 
1972  and  2  November  1969.  Their  data  show  proton  and  electron  fluxes 
measured  in  several  energy  bands  as  the  satellites  traversed  the  polar 
cap  at  altitudes  of  several  hundred  kilometers.  During  the  first  half¬ 
day  of  the  2  November  1969  event,  the  fluxes  of  1.7,  13.0,  and  24.5  MeV 
protons  exhibited  spatial  fluctuations  up  to  about  a  factor  of  3  within 
the  cap.  No  data  were  obtained  for  the  50  to  100  MeV  protons,  which 
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ionize  the  important  40  to  50  km  region.  Because  data  on  the  lower 
altitudes  are  insufficient  and  the  chemical  coefficients  are  imprecise, 
it  is  pointless  to  attempt  a  detailed  calculation  of  the  propagation 
uncertainty  due  to  the  factor-of-3  spatial  "jitter"  in  proton  flux. 

The  ion  density  varies  only  as  the  square  root  of  the  production  rate, 
which  is  proportional  to  the  incident  flux.  Thus,  the  spatial  in- 
homogeneitv  of  the  flux  will  probably  cause  no  more  than  a  factor  of 
2  inhomogenei ty  in  the  refractive  index.  As  discussed  below,  addi¬ 
tional  variations  will  occur  if  reaction-rate  coefficients  are  non- 
uniform  across  the  cap. 

Spatial  variations  in  proton  flux  seem  to  diminish  as  the  PCA 
develops.  Data  taken  on  3  November  1969  indicate  a  more  uniform  flux 
than  that  for  2  November.  Moreover,  data  given  by  Reagan  and  Watt 
[1976]  show  that  proton  fluxes  measured  at  several  energies  between 
1  and  100  MeV  were  quite  uniform  across  the  cap  on  4  August  1972 — the 
second  day  of  that  event. 

These  few  measurements  indicate  that  only  slight  lateral  in¬ 
homogeneities  in  the  refractive  index  are  caused  by  inhomogeneities 
in  proton  flux  once  a  PCA  has  had  a  day  to  develop.  Even  earlier, 
the  refractive  index  inhomogeneitv  due  to  flux  inhomogeneitv  is  prob- 
ablv  no  worse  than  a  factor  of  2.  Sensitivity  studies  show  that  no 
intolerable  uncertainties  in  predicted  ELF  signal  strengths  would 
occur  even  if  the  conductivity  height  profile  were  uncertain  by  a 
factor  of  2  over  the  entire  first  Fresnel  zone.  In  fact,  the  satel¬ 
lite  data  indicate  that  transverse  variations  in  proton  flux  have 
characteristic  scales  of  only  hundreds  of  kilometers,  and  would  be 
averaged  over  the  enormous  ELF  Fresnel  zones  (see  Figs.  5  and  6  below). 

Having  discussed  inhomogeneities  due  to  lateral  structure  in  the 
proton  flux — i.e.,  those  that  would  occur  even  if  reaction-rate  coef¬ 
ficients  were  uniform  across  the  cap — we  conclude  that,  in  the  two  PCAs 
examined,  such  inhomogeneities  apparently  would  not  have  seriously 
degraded  the  interpretation  of  a  transpolar  ELF  propagation  experi¬ 
ment.  However,  lateral  inhomogeneities  can  be  caused  bv  another 
mechanism.  The  chemical  coefficients  governing  ionization  caused  bv 
a  given  proton  flux  could  vary  across  the  cap.  If  such  a  variation 
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were  large,  lateral  gradients  in  the  refractive  index  could  result  even 
if  tlie  flux  were  spatially  uniform.  To  evaluate  the  implications  of 
such  chemistry-induced  inhomogeneities,  air  chemistry  research  must 
estimate  the  variation  of  key  reaction-rate  coefficients  across  the 
cap . 


EFFEGT  OF  POLAR  CAP  BOUNDARY  ON  ELF  PROPAGATION 

Since  transpolar  ELF  propagation  depends  on  the  size  of  the 
proton- irrad iated  region,  it  is  affected  by  polar  cap  boundary  vari¬ 
ability.  Satellite  data  on  the  geographic  dependence  of  the  proton 
flux  are  given  by  Imhoff  et  al.  [1976].  For  the  energetic  protons 
of  interest  here,  the  effective  boundary  of  the  cap  typically  occurs 
within  invariant  latitudes  of  65  deg  +  5  deg;  i.e.,  the  linear  dimen¬ 
sions  of  the  cap  vary  up  to  about  5  deg/(90  deg  -  65  deg)  ~  20  per¬ 
cent,  whereas  the  corresponding  area  varies  up  to  about  40  percent. 

Figures  5  and  6  show  the  first  Fresnel  zone  for  75  Hz  propagation 
from  WTF  to  Tromso  and  Thule,  respectively.  The  shaded  area  is  the 
region  irradiated  by  energetic  protons  at  1400  UT  on  4  August  1972 
[Imhoff  et  al.,  1976].  As  noted  above,  the  boundaries  of  the  polar 
cap  can  vary  by  a  few  degrees  in  latitude.  Here  we  investigate  the 
sensitivity  of  transpolar  propagation  to  the  location  of  the  effective 
boundaries . 

It  is  convenient  to  characterize  polar-cap-boundary  effects  as 
either  longitudinal  or  transverse.  Longitudinal  effects  depend  on 
the  length  of  the  great-circle  propagation  path  exposed  to  the  PGA, 
whereas  transverse  effects  depend  on  the  fraction  of  the  I'resnel-zone 
cross  section  filled  by  the  PCA.  For  a  laterally  uniform  cap  with 
sharp  boundaries,  longitudinal  effects  are  governed  by  the  intersec¬ 
tions  of  the  cap  boundary  with  the  great-circle  propagation  path 
(see  Figs.  5  and  6);  transverse  effects  are  governed  by  the  inter¬ 
sections  of  the  boundarv  with  the  line  that  normally  bisects  the  great- 
circle  path  and  terminates  at  the  widest  points  of  the  first  Fresnel 
zone.  The  length  of  this  line  is  2d,  where 


d  = 


(9) 
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Irradiated  polar  cap 


irssuaq 


is  the  maximum  half-width  of  the  Fresnel  zone,  x  is  the  great-circle 
pathlength,  X  is  the  wavelength,  and  the  propagation  constant  ~ 

1.15  to  1.5  for  75  Hz  propagation 

Longitudinal  boundary  effects  are  relatively  straightforward  to 
analyze.  For  example,  consider  the  paths  from  the  WI'F  to  Tromso  and 
Thule  in  Figs.  5  and  6.  We  assume — subject  to  later  deliberation — 
that  (1)  the  cap  fills  the  Fresnel-zone  cross  section,  so  transverse 
effects  are  minimal;  and  (2)  the  transmitter  is  outside  the  cap’s 
boundaries  and  the  receiver  is  within.  The  latter  assumption  is  valid 
for  the  WTF-to-Thule  path,  but  could  be  violated  for  WTF-t o-Tromso 
should  the  boundary  occur  at  a  latitude  higher  than  indicated  in  Fig.  5. 

Subject  to  these  assumptions,  uncertainties  in  the  boundary  loca¬ 
tion  can  affect  only  the  PCA-induced  attenuation  A  (but  not  the  excita¬ 
tion),  which  is  given  by 


A  =  (B  -  Bq)1.  dB  ,  (10) 

where  3  is  the  attenuation  rate  under  the  disturbed  cap,  3Q  is  the 
ambient  attenuation  rate,  and  I.  is  the  exposed  pathlength.  The 
fractional  uncertainty  in  attenuation  AA/A  caused  by  an  uncertainty 
in  exposed  pathlength  Ah,  due  to  an  imprecisely  known  boundary  loca¬ 
tion,  is  simply 


AA  _  AL 
A  ~  lT  • 


(ID 


In  short,  without  continuously  monitoring  the  extent  of  the  cap, 
the  boundary  location  can  be  estimated  to  within  +5  deg  of  latitude, 
i.e.,  to  within  about  +500  km.  1.  «  5  Mm  for  a  transpolar  path  such 
as  WTF  to  Tromso,  and  about  2.5  Mm  for  WTF  to  Thule.  The  correspond¬ 
ing  fractional  uncertainties  in  PCA-induced  attenuation  art*  +10  ami 
+20  percent,  respectively.  The  absolute  value  of  the  longitudinal 
effect  of  a  500  km  uncertainty  in  boundary  location  would  be  about 
+0.5  dB  uncertainty  in  calculated  field  strength  for  a  relatively 
strong  event  (1  dB/Mm  excess  attenuation).  These  numbers  could  double 
if  the  receiver  were  outside  the  cap,  because  boundary  errors  would 
occur  at  each  end  of  the  path. 
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A  more  important  error  could  result  if  the  receiver  were  so 
near  the  boundary  that  it  was  uncertain  whether  the  disturbed  or  am¬ 
bient  value  should  be  used  for  the  excitation  factor.  For  a  strong 
event,  the  difference  between  ambient  and  disturbed  excitation  at 
the  receiver  could  be  ~1 . 5  dB  [Field,  1969].  Thus ,  the  status  of 
the  receiver  site  must  be  established  either  by  sounding  the  iono¬ 
sphere  or  locating  the  receiver  so  that  no  reasonable  doubt  exists 
as  to  whether  it  is  inside  or  outside  the  irradiated  cap. 

Transverse  effects  are  more  complicated  than  longitudinal;  they 
can  be  analyzed  with  either  scattering  theory  or  the  integral  wave 
equation.  Field  and  Joiner  [1979]  used  the  latter  approach  to  obtain 
results  that,  although  not  derived  explicitly  for  a  polar  cap  model, 
clarify  the  problem  at  hand.  Specifically,  Field  considered  a  dis¬ 
turbance  having  the  transverse  dependence 

-y2/ (Ay) 2 

* 

wtiere  y  is  the  distance  from  the  direct  propagation  path  and  Ay  is 
a  measure  of  the  transverse  size  of  the  disturbance.  It  was  shown 
that  the  ratio  R  of  the  actual  field-strength  reduction  to  that  com¬ 
puted  by  assuming  the  disturbance  to  be  infinitely  wide  (i.e.,  by 
using  the  WKB  approximation  and  thus  ignoring  transverse  gradients) 
is  simply 


If  the  disturbance  is  wider  than  a  Fresnel  zone  (d/Ay  <•'  1),  then 
R  ->  1 ,  the  effects  of  finite  transverse  size  are  negligible,  and  the 
WKB  result  is  recovired.  For  Ay  ^  d,  R  is  smaller  than  unity,  indi¬ 
cating  that  PCA- induced  attenuation  computed  including  transverse 
effects  is  less  than  that  computed  by  simply  integrating  along  the 
direct  path  (e.g.,  as  in  F,q .  (10)). 


Figure  7  shows  R  as  a  function  of  Ay  and  Av/d  for  a  6  Mm  path 
and  frequencies  of  45  Hz  and  75  Hz.  Transverse  effects  would  he 
negligible — and  the  transverse  location  of  the  cap  boundaries  unim¬ 
portant —  if  |r|  were  unity  and  the  phase  of  R  were  zero.  Tf  the 
disturbance  fills  at  least  75  percent  of  the  Fresnel-zone  cross  sec¬ 
tion,  then  | R ]  >  0.9  and  phase  R  <  10  deg,  which  indicates  that  boundary 
effects  are  unimportant  if  Ay/d  >  0.75.  Figures  5  and  h  demonstrate 
that,  for  the  paths  shown,  more  than  80  percent  of  the  Fresnel-zone 
cross  section  is  filled.  The  error  made  by  ignoring  edge  effects  is 
about  5  percent  (Fig.  7).  Of  course,  transverse  boundary  effects 
could  be  much  greater  for  other  paths — e.g..  Ay/d  would  be  about  0.5 
for  WTF  to  Scotland.  In  such  cases,  the  location  of  the  cap  boundary 
must  be  known  to  correctly  calculate  the  signal.  The  simplest  pro¬ 
cedure  would  be  to  avoid  paths  passing  near  the  edge  of  the  cap  and 
ignore  transverse  boundary  effects. 


L 


I 


Ay/d 


0  500  1000  1500  2000  2500  3000 


Ay  (km) 

Note  d  ”  Fresnel  rone  half-width,  .iy  “  effective  half-width  of  disturbance 

Figure  7.  Ratio  of  actual  propagation  anomaly  to  that  computed 
ignoring  transverse  boundary  effects. 
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Phase  R  (deg) 


SECTION  4 


SICNAL  AND  NOISE  PROCESSING 


Here  we  establish  the  amount  of  signal  and  noise  processing  needed 
to  accurately  estimate  the  transpolar  WTF  signal.  The  tradeoff  is 
difficult  because  the  signal  integration  time  must  be  long  enough  to 
yield  a  good  postprocessing  s i gna I - t o-no i se  ratio  (SNR),  but  short 
relative  to  temporal  changes  in  a  PCA.  One  troublesome  aspect  is 
that  strong  PCAs--which  best  simulate  nuclear  environments — cause 
the  lowest  SNRs  and  thus  require  tlu  longest  integrations.  After 
estimating  atmospheric  noise  and  the  expected  ambient  and  disturbed 
signal  strengths  at  two  possible  receiver  sites,  wo  determine  the 
integration  times  needed  to  resolve  the  PCA- induced  change  in  signal 
amp  1  i t  ude . 


EXPECTED  SIGNAE  AND  NOISE 

We  assume  that  the  WTF  signal  will  he  received  at  Thule,  Greenland 
am!  Trotnso,  Norway,  whose  positions  are  shown  in  Figs.  r>  and  b.  We 
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Table  i .  Assumed  waveguide  parameters 
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Table  4.  Calculated  signal  amplitudes. 


Normal 

Disturbed 

Location 

(dBH) 

(dBH) 

Tromsci 

-156 

-164 

Thule 

-150 

-155 

NOTE:  Values  are  in  decibels 

relative  to  1  A/m. 


Keans  and  Griffiths  [1974]  as  well  as  Davis  and  Mvers  [197b] 
report  measurements  of  effective  atmospheric  noise  in  Norway  after 
excision  or  clipping  of  the  larger  impulses.  The  former  cite  a 
median  of  -14?  dlill,  with  lows  of  -148  dBH  and  highs  of  -138  dBH. 

They  observed  a  diurnal  variation  and  found  a  2  dB  decrease  in  noise 
in  the  carle  morning  in  Norway.  Davis  and  Myers  [197b],  who  investi¬ 
gated  impulse  clipping  at  various  levels,  record  slightly  higher 
<  fleet i ve  noise  levels,  with  a  median  of  nearly  -141  dBH  during  the 
.summer  in  Norway.  Thee  estimate  that  processing  reduced  effective 
noise  be  over  10  dB.  Without  noise  processing,  therefore,  their  mea¬ 
sured  level  would  have  been  about  -130  dBH. 

We  compromise  between  the  above  measurements,  using  -142  dBH  for 
the  expected  atmospheric  noise  level  at  tile  receiver  it  clipping  is 
Used.  We  estimate  noise  to  he  10  d B  higher  in  the  absence  of  clipping. 
Table  i  shows  expiated  signa 1 -t o-ito iso  ratios  obtained  he  combining 
these  noise  levels  with  tile  signal  strengths  shown  in  Table  4.  We 
havi  used  tiie  same  noise  levels  for  ambient  and  disturbed  condit  ions 
at  Ironso,  which  is  mar  tin  edge  ot  t  lie  polar  cap  and  so  receives 
a  t  mo  spiie  r  i  c  no  i  e  via  nonpolar  propagation  paths  from  most  thunder- 
•  term  centers.  I'iiule,  on  the  other  hand,  is  well  within  the  cap,  and 
we  have  assumed  that  a  i’CA  would  suppress  signal  and  noise  about 
eipiai  l  V.  Unis ,  tin  disturbed  and  undisturbed  SNKs  are  taken  to  tie 
t  he  same  t  here  . 
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Table  5.  Expected  SNRs  under  normal  and  disturbed  conditions. 


Normal  (dB) 

Disturbed 

(dB) 

Locat ion 

Clipped  Unclipped 

Noise  Noise 

Cl  i pped 

Noi  se 

Une 1 i pped 
No  i  se 

Tromso 

-14  -24 

-22 

-32 

Thule 

1 

oc 

1 

M 1 

CO 

-8 

-18 

NOTE:  Values  pertain  to  a  1  Hz  bandwidth. 


INTEGRATION  OF  SIGNAL 

For  verifying  the  predictive  codes,  the  signal  must  be  measured 
precisely  enough  to  resolve  the  change  caused  by  a  PGA.  For  a  very 
strong  event.  Table  4  shows  that  an  8  dB  decrease  in  tiie  signal  would 
be  expected  at  Tromso.  In  such  a  case,  the  signal  would  have  to  be 
measured  to  within  about  2  dB  before  and  during  the  event.  The  mea¬ 
surement  during  the  event  is  the  more  difficult,  beet,  se  of  the  re¬ 
duced  SNR.  Events  small  enough  to  cause  a  signal  reduction  of  just 
a  few  decibels  would  require  measurement  precision  of  1  dB  or  less. 
However,  that  necessity  is  partially  mitigated  by  the  greater  SNR 
during  a  moderate  than  a  severe  event. 

Table  5  shows  SNR  in  a  1  Hz  bandwidth,  which  is  equivalent'  to 
a  1  sec  integration  time.  Since  these  ratios  are  very  low,  integra¬ 
tion  times  much  longer  than  1  sec  are  needed  to  measure  the  signal 
with  acceptable  statistical  stability.  Two  kinds  of  integration  can 
be  considered — coherent  and  noncoherent.  Coherent  integration  re¬ 
quires  that  signal  phase  be  constant  over  the  period  of  integration, 
or,  alternatively,  that  estimates  of  phase  be  available  so  variations 

can  be  compensated.  The  usual  practice  is  to  integrate  coherent lv  as 

long  as  possible,  subject  to  both  conditions  regarding  phase,  then  in¬ 
tegrate  noncoherent ly  for  further  improvement  in  postprocessing  SNR. 
During  the  period  of  coherent  integration,  the  SNR  increases  linearly 
with  time.  Noncoherent  integration  improves  SNR  at  a  lesser  rate — 
a  rule  of  thumb  being  that  the  improvement  is  proportional  to  the 

square  root  of  the  integration  time.  Here  we  estimate  the  expected 

SNR  somewhat  more  accurately. 


Wo  first  assume  that  the  signal  is  integrated  coherently  and 
compute  the  likelihood  that  noise  causes  the  measured  signal  to  exceed 
tile  expected  signal  by  various  amounts.  Specifically,  we  compute  the 
power  level  that  would  be  exceeded  by  chance  on  no  more  than  10  per¬ 
cent  of  independent  attempts  to  measure  the  signal  under  specific 
conditions . 

That  likelihood  can  be  expressed  in  terms  of  the  noncentral  chi- 
square  distribution  with  two  degrees  of  freedom.  Let  the  expected 


complex  signal  amplitude  be  rO  e*^  = 
uniformly  distributed  in  phase  <p  with  standard  deviation  a.  The  like 


m  +  m  ,  where  the  noise  is 
x  y’ 


lihood  that  a  measurement  of  the  power  would  be  equal  to  or  less  than 
2  2 

R  O  is  given  by  the  integral 


//< 


(2ttc2)  ^  exp 


(x-m)^+(y-m)^_l 


2a 


dx  dy  =  P(R2 | 2,  r2l  , 


(15) 


where  A  denotes  the  area  of  a  circle1  of  radius  RO  centered  at  the 

2 1  2 

origin,  and  P(R  |2,  r  )  is  the  cumulative  distribution  function  of 

the  noncentral  chi-square  distribution  with  two  degrees  of  freedom 

,  2  2 
and  noncentrality  parameter  r  .  Tn  Eq .  (15),  we  may  identify  r  as 

? 

the  expected  SNR  in  a  l  II z  bandwidth  and  R“  as  the  corresponding 
power  realized  in  a  measurement. 

If  a  is  the  rms  noise  amplitude  in  a  1  Hz  bandwidth,  the  rms 
noise  amplitude  resulting  from  T  sec  of  coherent  integration  would 
be  Ov^T,  whereas  the  integrated  signal  amplitude  would  be 

rOT  e ' ^  =  m  T  +  im  T  . 

x  v 

The  likelihood  that  the  integrated  measurement  would  indicate  a  power 
2  2  2 

no  greater  titan  R^O  T*'  is,  by  analogy  with  Eq .  (15),  given  by 


//<  2”2T)“! 


exp 


(x  -  ™XT) * 


)2  +  (y  -  m^T)2" 
2o2T 


dx  dy  =  P(R2T|2,  r2T)  , 


(16) 


where  A  is  now  the  area  within  a  circle  of  radius  RaT. 

2  . 

Tables  of  the  noncentral  chi-square  distribution  P(x  |v,  A)  are 
2 

available.  For  x  and  A  large  (>3),  however,  the  cumulative  distribu¬ 
tion  function  can  be  approximated  with  the  normal  distribution  function 


P(x2|v,  A) 


(17) 


where 


,  2  1/3  2  / 1  +  b  \ 

(x  /e) _ -L1  - 

_ _ * 

yi  a  +  b)/g 


(18) 


and  where  g  =  V  +  A  and  b  =  A/g  [Abramowitz  and  Stegun,  1972].  Equa¬ 
tions  (16)  through  (18)  were  used  to  compute  Fig.  8,  which  shows, 

2  i  2  2  2 

for  P(R  T | 2 ,  r  T)  =  0.9,  how  R  approaches  r  with  increasing  coherent 
integration  time  T.  The  SNR  ratio  in  a  1  Hz  bandwidth  is  shown 
parametrically  on  the  curves. 

Figure  8  shows  envelopes  of  90  percent  confidence  that  the  ratio 
R/ r  of  the  measured  signal  to  the  mean  signal  will  not  be  greater 
than  shown  on  the  ordinate.  A  set  of  similar  curves  would  show  the 
envelopes  of  90  percent  confidence  that  the  inverse  of  this  ratio  will 
not  he  greater  than  shown  on  the  ordinate.  The  curves  on  Fig.  8  show 
envelopes  of  80  percent  confidence  that  the  quantity  | 20  log  R/r]  will 
thus  not  he  greater  than  the  amount  indicated. 

To  illustrate  the  use  of  Fig.  8,  consider  normal  conditions  at 
Tromso  with  noise  clipping  (SNR  =  -14  dP,  from  Table  5).  The  figure 
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percent  confidence  levels  versus  coherent  int 


shows  that  1  hr  of  coherent  integration  would  give  80  percent  con¬ 
fidence  that  the  measured  signal  is  within  +1  dB  of  the  mean  signal. 
For  highly  disturbed  conditions  at  Tromso  with  noise  clipping 
(SNR  =  -22  dB  from  Table  5),  Fig.  8  siiows  that  2  x  10^  sec  (about 
5.5  hr)  of  coherent  integration  is  needed  to  achieve  ±1  dB  accuracy 
witli  80  percent  confidence.  The  importance  of  noise  clipping  is 
demonstrated  by  using  SNRs  for  unclipped  noise  from  Table  5  in  con¬ 
junction  witii  Fig.  8.  For  normal  conditions  at  Tromso  (SNR  =  -24  dB), 
nearly  10  hr  of  coherent  integration  is  needed  for  80  percent  confi¬ 
dence  of  ±1  dB  accuracy.  For  disturbed  conditions  (SNR  =  -32  dB), 
the  requirements  are  so  severe  that  the  80  percent  confidence  envelope 
is  off-scale  on  Fig.  8.  Thus  noise  clipping  is  essential. 

This  theory  is  valid  either  if  the  signal  phase  is  constant,  or 
if  its  variation  is  known  and  can  be  compensated.  Less  efficient  non¬ 
coherent  integration  must  be  used  if  neither  condition  is  met.  In 
such  a  case,  coherent  integration  would  be  employed  as  long  as  pos¬ 
sible — e.g.,  as  long  as  phase  could  be  assumed  to  vary  by  less  than, 
say,  1  rad;  or  as  long  as  permitted  by  the  particular  receiver  used. 

At  the  end  of  each  such  interval,  the  coherently  integrated  power 
would  be  measured;  and  the  results  of  some  number  N  of  such  measure¬ 
ments  would  be  averaged  to  obtain  a  single  estimate  of  the  received 
power . 

If  the  complex  amplitude  of  the  ktli  measurement  of  the  coherently 
integrated  signal  (plus  noise)  is  +  iY^,  the  result  of  summing  the 
power  for  N  such  measurements  is 


2 

x 


2  2 

<xk  +  V 


k=l 


Tins  summation  again  tins  a  noncentral  chi-square  distribution,  but 

now  witii  2N  degrees  of  freedom.  Tiie  likelihood  that  the  averaged 

2  2,2 

power  measurements  do  not  exceed  R  is  thus  given  by  P(R  T(,  1 2N,  r  T  ), 

where  T  is  the  duration  of  the  coherent  integration  and  N  is  the 
c 

number  of  coherent  measurements;  i.e.,  N  =  T/T^ . 
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Figure  9,  analogous  to  Fig.  8,  shows  the  80  percent  confidence 

envelopes  for  noncoherent  averaging  of  measurements  made  with  hour- 

•k 

long  coherent  integrations.  Note  that ,  for  a  given  integration  time, 
the  signal  uncertainty  indicated  in  Fig.  9  is  somewhat  greater  than 
that  shown  in  Fig.  8  for  coherent  integration. 

Figure  10  presents  the  information  from  Fig.  9  so  as  to  illustrate 
the  effect  of  integration  time  on  resolution  of  changes  in  received 
power.  The  horizontal  lines  correspond  to  mean  SNR  in  a  1  Hz  band¬ 
width,  whereas  the  shaded  regions  surrounding  eacli  line  give  the  80 
percent  confidence  envelopes.  The  labels  on  the  horizontal  lines 
correspond  to  the  values  for  clipped  noise  given  in  Table  5. 

Consider,  for  example,  a  receiver  at  Tromso  where  the  mean  ambient 
SNR  (in  1  Hz)  is  -14  dB.  A  very  strong  PCA  would  change  it  to  -22  dB — 
an  8  dB  separation  between  the  disturbed  and  undisturbed  values. 
Statistically  valid  resolution  of  that  change  requires  an  integration 
time  long  enough  that  the  80  percent  confidence  values  corresponding 
to  the  two  mean  signal  levels  are  well  separated.  The  graph  shows 

4 

that  3  hr  ( — 10  sec)  of  integration  would  resolve  the  normal  signal 
to  within  +0.6  dB  witli  80  percent  confidence,  and  the  disturbed  signal 
to  within  +1.6  dB.  The  resulting  ~2  dB  uncertainty  in  the  difference 
between  the  disturbed  and  undisturbed  signals  is  probably  tolerable 
for  an  "8  dB"  event . 

Much  longer  integration  times  are  needed  to  resolve  smaller  events. 

Consider,  for  example,  a  moderate  "4  dB"  PCA  that  reduces  the  mean 

4 

Tromso  SNR  from  -14  to  -18  dB.  A  10  sec  integration  would  result  in 
a  +1  dB  uncertainty  in  disturbed  signal.  Combined  with  the  +0.6  dB 
uncertainty  in  normal  signal,  that  value  would  yield  a  total  uncertainty 
of  ~1 . 6  dB — 40  percent  as  large  as  the  4  dB  change  being  monitored, 
and  probablv  too  large  for  code  verification.  A  total  uncertainty 
less  than,  say,  1  dB  would  be  acceptable  in  this  case. 

Figure  10  shows  that  an  integration  time  of  nearly  10  hr  would 
be  needed  to  achieve  +0.6  dB  and  +0.2  dB  uncertainty  on  the  disturbed 

* 

This  choice  of  coherent  integration  time  is  compatible  with  t he 
capabilities  of  receivers  considered  for  the  experiment. 
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Figure  9.  80  percent  confidence  levels  versus  time:  noncoherent  integration 

of  hour-long  coherently  integrated  samples. 


percent  confidence  envelopes  versus  time:  noncoherent 
egration  of  hour-long  coherently  integrated  samples. 


and  normal  signals,  respectively.  Such  long  times  are  troublesome 
because  the  ionosphere — and,  hence,  the  signal  level — could  well 
change  during  the  measurement.  Figure  9  shows  that  some  reduction 
in  the  integration  time  would  be  possible  if  coherent  integration 
could  be  used  throughout  the  entire  measurement.  Otherwise,  a  rela¬ 
tively  strong  event  is  needed  if  a  statistically  reliable  measurement 
of  the  change  in  signal  level  is  to  be  obtained. 
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